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Abstract 

We develop a framework to track the structure of temporal networks with a signal processing approach. The 
method is based on the duality between networks and signals using a multidimensional scaling technique. This 
enables a study of the network structure using frequency patterns of the corresponding signals. An extension is 
proposed for temporal networks, thereby enabling a tracking of the network structure over time. A method to 
automatically extract the most significant frequency patterns and their activation coefficients over time is then 
introduced, using nonnegative matrix factorization of the temporal spectra. The framework, inspired by audio 
decomposition, allows transforming back these frequency patterns into networks, to highlight the evolution of 
the underlying structure of the network over time. The effectiveness of the method is first evidenced on a toy 
example, prior being used to study a temporal network of face-to-face contacts. The extraction of sub-networks 
highlights significant structures decomposed on time intervals. 


1 Introduction 

Many complex systems, whether physical, biological or social, can be naturally represented as networks, i.e., a set 
of relationships between entities. Network science [28] has been widely developed to study such objects, generally 
supported by a graph structure, by providing powerful tools, such as the detection of communities [12], in order 
to understand the underlying properties of these systems. Recently, connections between signal processing and 
network theory have emerged: The field of signal processing over networks mills] has been introduced with the 
objective of transposing concepts developed in classical signal processing, such as Fourier transform or wavelets, 
in the graph domain. These works have led to significant results, among them filtering of signals defined over a 
network [31] or multiscale community mining using graph wavelets [33] . These approaches benefit from the natural 
representation of networks by graphs, enabling the use of the comprehensive mathematical understanding of graph 
theory. Another approach has been also considered, defining a duality between graphs and signals: methods have 
been developed to transform graphs into signals and conversely, in order to take advantage of both signal processing 
and graph theory. Hence, mapping a graph into time series has been performed using random walks [HIEIIIS] 
or deterministic methods based on classical multidimensional scaling [23] [30] . This latter approach has been the 
topic of several extensions in [22], in order to build a comprehensive framework to link frequency patterns of the 
so-obtained signals with network structures. 

Studies mainly focused on the analysis of static networks, potentially aggregated over time, but without time 
evolution. However, the considered systems are most of the time not frozen: vertices and edges appear and disappear 
over the course of time. Aggregating the networks over a time interval gives insight of the underlying mechanisms, 
but often does not provide the actual dynamic sequence of appearance and disappearance of edges: two edges, active 
one after the other, will be considered simultaneous in the temporally-aggregated network. Given the importance 
of knowing such dynamics, for instance in topics such as epidemic spread or communication networks, and thanks 
to the recent availability of many data sets, a temporal network theory has recently appeared piiij. enabling a 
deeper understanding of underlying mechanisms. Several studies proposed an extension of the methods developed 
for static networks to the temporal case: we can cite for instance works on network modeling laiaiisiEi], detection 
of communities isaisaiiiiiiHiis], detection of temporal motifs [25] , visualization [36] , or more generally data mining 
of time-evolving sequences [38] . 

We propose in this article a new approach based on the duality between networks and signals to study temporal 
networks, introduced in [22]. The objective here is to follow the global structure of the network over time. One 
first contribution consists of an extension of the existing work to the temporal case, which is naturally performed 
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by considering each time instant as a static network. This approach enables us to visually track temporal networks 
by following the frequency patterns associated to specific structures. The second contribution consists of using 
nonnegative matrix factorization (NMF) [26] to automatically extract and track the significant frequency patterns 
over time. Details about the reconstruction how components can be transformed back into networks are also 
provided. 

Preliminary versions of this work have been presented in different places: the principle of extension of the 
method to the temporal case is introduced in m, as well as the visual tracking of frequency patterns representing 
structures. In min], this approach has been used to study a bike sharing system in Lyon, which is not mentioned 
in this article. Finally, the idea of the decomposition using NMF has been suggested in [20]. This paper extends 
however those previous works by detailing a comprehensive framework and setting out consistent arguments for the 
validity of the method in the context of analysis of temporal networks. 

The paper is organized as follows. Sectionj^recalls the duality between static networks and signals, as introduced 
in m and extended in [22]. Section]^ gives an extension of the transformation to the temporal case, and introduces 
a toy temporal network highlighting specific network structures. Section [^ shows how spectral analysis is used to 
track the structure of the temporal network over time, while Section describes a decomposition of a temporal 
network into sub-networks by applying nonnegative matrix factorization to the matrix of spectra over time. For 
both sections, an illustration with the toy temporal network is provided. Finally, Section [^ applies the developed 
framework to a real-world temporal network from the SocioPatterns project, describing face-to-face contact between 
children in a primary school. 

Notations Matrices are denoted by boldface capital letters, their columns by boldface lowercase letters, and 
their elements by lowercase letters: for M G M = [mi,...,mB] = Tensors are 

represented by calligraphic capital letters: for T G T = Operators are represented by 

calligraphic boldface capital letters: 

2 Duality between static networks and signals 

2.1 Transformation from static networks to signals 

Duality between networks and signals has been introduced to analyze networks structures using signal processing 
tools. Shimada et al. m proposed a method to transform static networks into a collection of signals using mul¬ 
tidimensional scaling. We proposed in [22] a comprehensive framework to transform graphs into a collection of 
signals, based on this work. We recall in this section the main principles of this framework, which will be used 
in the following to study temporal networks. We consider in the following networks described by unweighted and 
undirected graph with N vertices. 

In [30], Classical MultiDimensional Scaling (CMDS) [T] is used to transform a graph into signals by projecting 
the N vertices of the graph in a Euclidean space, such that distances between these points correspond to relations 
in the graph. The graph is fully described by its adjacency matrix A G whose elements aij are equal to 1 if 

vertices i and j are linked, and 0 otherwise. The transformation consists in applying CMDS to a matrix encoding 
distance between vertices of a graph, noted A = and defined for two vertices i^j of the network by: 

0 if i = jf 

1 if aij = 1 and i j (1) 

w > 1 if aij = 0 and i j 

As discussed in [22], we choose w = 1 This definition focuses on the presence (denoted by a distance equal 

to 1) or the absence (denoted by a distance equal to w) of an edge between two vertices. Hence, the distance of 
two vertices in the graph, often defined as the length of the shortest path between the two vertices, has no direct 
impact on the matrix A: two pairs of unlinked vertices will have a distance equal to w, whether they are close or 
not in the graph. Applying CMDS on the distance matrix A leads to a collection of points, corresponding to the 
vertices, in a Euclidean space The considered signals, or components, correspond to the coordinates, for 

each dimension, of the points. The collection of signals is denoted by X G where C is the total number of 

components, and thus is equal to N — 1 for the full representation. The columns Xc represent the c-th signal, with 
c G {1,... C}, and are indexed by the vertices. The transformation is denoted by 'T : T[A] = X. 
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2.2 Inverse Transformation 

Transforming back signals to a graph has to take into account the nature of the signals, as they are a representation 
of a specific network. The inverse transformation must hence preserve the original topology of the underlying graph. 
By construction of the collection of signals X, the perfect retrieval of the underlying graph is easily reachable, by 
considering the distances between each pair of point: As built using CMDS, these distances represent the distance 
matrix A, and the adjacency matrix of the graph is directly obtained. However, when X is degraded or modified, 
e.g. by taking C < X — 1, the distances are no longer directly the ones computed between vertices, even if they 
stay in the neighborhood of these distances. We proposed in [22] to take into account the energy of components 
to improve the reconstruction, as well as prior information about the original graph. If the distance between two 
vertices i and j in a high-energy component is high, it means that the two vertices are likely to be distant in the 
graph. Conversely, if the distance in a high-energy component is low, then the two vertices are likely to be connected 
in the graph. 

Let X be a degraded collection of signals. The energies are normalized according to the energies of the original 
components, by multiplying the components Xc by the normalization factor Nc'. 


Nr = 


2 

nc 


W x' 
Z^n=l 

^ Z^n=l '^nc 

Then the distances are computed by using the energies as follows: 
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with a > 0, where Uc is the energy of component c, computed as Uc = and normalized such that 

Iklla = E^=l = C', with u the vector of energies for all components. The parameter a controls the importance 
of the weighting: if a is high, the high-energy components have a higher importance in the computation of distances 
compared to the low-energy components. Conversely, if a is small, the importance of high-energy components is 
diminished. In particular, a = 0 gives the standard reconstruction. The distances are then thresholded, by retaining 
the edges corresponding to the smallest distances. The threshold is chosen in order to recover the same amount of 
edges than in the original network. This inverse transformation is denoted 7 1 ^[X] = where is 

the adjacency matrix of the reconstructed graph. 


2.3 Indexation of signals 

As the signals are indexed by the vertices, the order in which we consider them in the transformation is essential 
to study some aspects of the signals, especially when using spectral analysis of the signals. Ordering randomly the 
vertices does not change the value assigned to each vertex, but would lead to abrupt variations in the representation 
of signals: Specific frequency properties, clearly observable in signals, will no longer be visible. Unfortunately, the 
suitable ordering is usually not available, especially when dealing with real-world graphs. To address this issue, 
we proposed in [21] to find a vertex ordering that reflects the topology of the underlying graph, based on the 
following assumption: if two vertices are close in the graph (by considering for instance the length of the shortest 
path between them), they have to be also close in the ordering. Details of the algorithm and results about the 
consistency between the obtained vertex ordering and the topology of the graph are covered in ED. 

2.4 Spectral analysis of signals 

Spectral analysis is performed using standard signal processing methods: Let a collection X of C signals indexed by 
N vertices. The spectra S G give the complex Fourier coefficients, whose elements are obtained by applying 

the Fourier transform on each of the C components of X: 

Sc = ^[Xc] (4) 

estimated, for positive frequencies, on F = X i bins, being the Fourier transform and cG {1,...,C}. 

From the spectrum S', the following features are obtained for each frequency of each component: 

• the magnitudes AT, which read as rricf = |sc/| 
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Table 1: Generation of the Toy Temporal Network: Probabilities to have an edge at time t 



eGEP 

e^EP 

e G Et-i 

0.99 

0.8 

e ^ Et-1 

0.2 

0.01 


• the energies which read as Ccf = \scf\‘^ 

• the phases which read as (/)c/ = arg(sc/) 

The matrices M and E are studied as a frequency-component map, exhibiting patterns in direct relation with 
the topology of the underlying graph. The phases of signals ^ are used in the inverse Fourier transformation, when 
the collection of signals has to be retrieved from M or E. 


2.5 Illustrations 

Figure shows illustrations of the transformation of a graph into a collection of signals, on several instances of 
network models. These illustrations show the connection between graph structure and the resulting signals after 
transformation. Regular structures (Figures p!a| and lb) highlight harmonic oscillations as components, while an 
organization in communities (Figure px) presents high-energy frequencies on the first components. Combination of 
both structures in the graph (Figure]^) is preserved in the frequency pattern. Finally, a random graph (Figure 
leads to a similar absence of structure in the corresponding frequency pattern. 


3 Extension to temporal networks 

3.1 Transformation of temporal networks into signals and back 

The description of temporal networks considered in this work consist in a discrete-time sequence of graph snapshots. 
The collection of all snapshots at the different times tG{0,...,T — 1}, with T the total number of time steps, can 
be represented by a graph adjacency tensor denoted A G We study here temporal networks where the 

edges are changing over time, keeping the same given set of vertices (possibly isolated). 

The extension of the method described in Section is directly achieved by applying at each time step the 
transformation on the corresponding static representation of the temporal network. We denote X G the 

collection of signals obtained from A. For each time step t, we have 

j^(t) ^ (5) 

As the number of vertices in the graph does not evolve, the number of components C, and then the number of 
frequencies F, is constant over time. In the case where the set of vertices evolves, the tensor is built by fixing 
the number of components as the maximal number of components over time, and by zero-padding the missing 
components. As for the indexation of vertices, the algorithm relabeling vertices according to the structure of the 
network is performed at each time step, leading to a different labeling of the vertices over timeQ 

Conversely, the inverse transformation is performed likewise, by applying the static inverse transformation on 
the collection of signals at time t: 


= T“i[X(‘)] (6) 

where represents the adjacency matrix of the reconstructed temporal network at time t. 

3.2 Illustration on a Toy Temporal Network (TTN) 

A toy temporal network is introduced for the sake of illustration. The temporal network consists of smooth 
transitions between different network structures. Starting from N unconnected vertices, the algorithm adds or 
removes edges at each time instant according to a probability depending both on the presence or not of the edge 

^The objective of this work here is to track how the structure of the temporal network evolves, regardless the labels of the vertices. 
The modification of the labeling indicates how the vertex evolves over time in the global structure of the network. This aspect is not 
discussed in the following, the labeling over time is kept only for purposes of reconstruction. 
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(a) 4-ring lattice 




(c) Stochastic block model with 3 communities, pw = 0.8 and pi) = 0.05 



(d) 4-ring lattice with 3 communities 



(e) Erdds-Renyi model with p = 0.4 


Figure 1: Illustrations on several instances of network models of the transformation of a network into a collection of 
signals. All networks have N = 100 vertices. (Left) Two-dimensional representation of the network. (Middle) First 
three highest energy components, and the arbitrarily chosen component 70, a low-energy one. (Right) Component- 
frequency map obtained after spectral analysis. The color represents the intensity, coded from white to black. The 
contrast is emphasized to expose the patterns. 
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(a) t = 19 (b) t = 39 (c) t = 59 



Figure 2: Adjacency matrix At of the Toy Temporal Network at different time instant. 


at the previous time, and on a prescribed network structure. Let Ef be the set of edges at time t, the set of 
edges of a prescribed network structure, and e = (i^j) the edge between the vertices i and j. Table gives the 
probability to have e in St according to the presence of e at the previous time (i.e. e G Et-i) and the presence 
of e in the prescribed network structure (e G E^). These probabilities are set up in order to generate a smooth 
transition between the initial structure and the prescribed network structure, and such that the structure at the 
end of the time interval is close to the one of the prescribed graph. 

Four networks are successively chosen as prescribed network structures, among those introduced in Section 
each structure being active for a time interval of 20 time steps: 


1. Random linkage of vertices (see Figure 


2 . 

3. 

4. 


Network with 3 communities (see Figure 
4-ring lattice (see Figure p!a[) 


4-ring lattice with 3 communities (see Figure Id) 


Figure shows the adjacency matrix of the temporal network at the end of each time interval. The prescribed 
structure is clearly visible thanks to a proper labeling of vertices: the blocks describe the communities, while the 
strong diagonal describes the regular lattice. These structures are nevertheless not completely defined, as noise due 
to the probabilities set in Table remains present. It allows to assess the robustness of the further analyses to 
noise. Figure shows the adjacency matrix of the temporal network aggregated over the 80 time steps: the four 
structures also appear using this representation. 

Figure plots three descriptors of the network at each time step, which would be a manner to analyze the tem¬ 
poral network using network theory. The number of edges (Figure [3a|), the average clustering coefficient (Figure [31^ 


6 

















































Avg. len. short, path Avg. clustering coefF. # of edges 



(a) Number of edges 




Figure 3: Descriptors of the Toy Temporal Network over time. The alternating shaded regions correspond to the 
four different periods, whose the label is given at the top of each region. 
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and the average length of shortest paths (Figure 3c) allow for an identification of the four time intervals: The 
number of edges increases in the period 2, to form the random structure, characterized by a low average clustering 
coefficient and low average length of shortest paths. Adding communities in period 2 increases the average cluster¬ 
ing coefficient, as expected. On the contrary, the regular structure in period 3 decreases the clustering coefficient 
and increases the average length of shortest paths, but much less in proportion than the decrease of the number 
of edges would affect. As expected, adding communities in this regular structure, as it is done in period 4, signifi¬ 
cantly increases the average clustering coefficient and slightly decreases the average length of shortest paths. These 
network-based descriptors give good intuitions on the underlying structure of the temporal network, but turn out 
to be inefficient to exactly characterize the structure. Furthermore, mixture of different structures, as it appears in 
this model, are not explicitly revealed. 


4 Spectral Analysis of Temporal Networks 

4.1 Temporal Spectra 

As previously, the extension of the spectral analysis of signals representing networks, defined in Section |T4| is simply 
achieved by considering independently each time step. We denote by 5 € the spectral tensor, where 

corresponds to the spectra obtained at time t. The corresponding matrices A4 and £ are respectively the temporal 
magnitudes and energies. 

As described in Sectionj^ the spectra are closely related to the network structures. In particular, the importance 
of high-energy components as well as low frequencies have been highlighted, for instance for the structure in 
communities. Looking at the marginals of the temporal energies or magnitudes over frequencies and components is 
hence expected to give hints about the evolution of the structure of the temporal network over time. In the following, 
we will focus on the marginal of the energies over the components, denoted Ec{f), and over the frequencies, denoted 
Ef{c). 

We also propose to use the spectral analysis to compare two network structures by computing the correlation 
between their spectra at each time step. The idea is to find among a set of parametric graph models the one that 
best fit the network at each time step. If we denote by Sm the spectra obtained after transformation of an instance 
Qrn of a prescribed graph model, we can compute a correlation coefficient between Sm and Generating 

several instances of the graph model gives us an average value of the correlation coefficient over several repetitions. 
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(a) Comparison with a network with k communities, for different values of k. The temporal network is highly correlated with a 
graph organized in communities during the period 2 and 4, as expected. 



(b) ^-regular lattice, for different value of k. The temporal network has the highest correlation with a 4-regular lattice during the 
period 3. 


Figure 5: Correlation between the temporal spectra at each time step and the spectra of networks with two specific 
network structures. The correlation is averaged over 20 repetitions. 

4.2 Illustration on the Toy Temporal Network (TTN) 

Figure shows Ec{f) and Ef{c) for respectively / G {1,2,3,35} and c G {0,1,2,70}. These two figures reveal 
the predominance of the first components and low frequencies, to track an organization in communities of the 
network: the low frequencies for the first components have a greater energies than for other types of structure, as 
already remarked in Section |T5l The tracking of the other types of structures is nevertheless not visible using this 
representation. 

The correlation between the temporal spectra of the toy temporal network and two network structures is studied. 
First, a structure in communities is observed, using a network model generating a random graph with a fixed number 
of communities. The comparison is done for a number of communities from 2 to 6, with 20 repetitions for each 
number of communities. Figure shows the average correlation: the correlation is maximal during the periods 2 
and 4, where the network is effectively structured in communities. During the period 2, the correlation is maximal 
when the network is structured in 3 communities, as expected. During the period 4, the number of communities is 
not clearly revealed, as the communities are not the only component of the network topology. 

In Figure |5b[ the temporal network is compared using the same method with a /^-regular lattice, for k equals 
to 2, 4, 6 and 8. If the correlations are lower than previously, we can nevertheless notice that in period 3, the 
temporal network is correlated with a 4-regular lattice, which is the structure set in the prescribed network during 
this period. 

The study of temporal spectra give hence insights about the structure of the underlying temporal network, but 
this approach is limited by the lack of knowledge, in the general case, of the parametric graph models composing 
the temporal network. To complement our analysis, we propose in the following to automatically find structures 
and their corresponding time activation periods that best characterize the temporal network. 
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5 Extraction of Frequency Patterns using Nonnegative Matrix Factor¬ 
ization 

5.1 Nonnegative Matrix Factorization (NMF) 

Nonnegative matrix factorization (NMF) [26] is a linear regression technique, used to decompose a nonnegative 
matrix V of dimension C x i.e. a matrix whose terms are greater or equal to zero, into the product of two 
nonnegative matrices W and H. NMF leads to a reduction of the dimensionality of data, by extracting in the 
columns of W patterns characterizing the data, and in the rows of H the activation coefficients of each pattern along 
the time. The number of extracted patterns is denoted K. A common approach to achieve such a decomposition 
consists of solving an optimization problem: 

= argmmD(V|TVff) (7) 

where is a dissimilarity measure. Fevotte et al. m proposed an algorithm to find a solution of the NMF where 
D is the /3-divergence, a parametric function with a simple parameter P which encompasses the Euclidean distance 
= 2 ), the generalized Kullback-Leibler divergence (/3 = 1) and the Itakura-Saito divergence (/3 = 0) as special 
cases. 

Regularization of the activation coefficients can be added in order to smooth them, with the assumption that 
there is no abrupt variations in the structure from one time step to the next one. The optimization problem is then 
defined as the minimization of the fitting term, defined in Eq. plus a term of temporal regularization: 

K T 

P{H) = Y^Y.^{hk(t-i)\hkt) ( 8 ) 

k=lt=2 

leading to: 


(FF*, iT*) = arg min D{V\WH) + -fP{H) (9) 

W H 

where 7 controls the regularization and is empirically fixed such that the activation coefficients highlight smoothness. 
In [9] and [8], smooth NME has been introduced for = 1 and f3 = 0. 

5.2 NMF on Spectra of Graphs 

Several approaches have been proposed to adapt NME to networks, either static [39] or temporal [T4|. In the latter 
approach, the adjacency matrix is represented as a tensor and is decomposed using nonnegative tensor factorization 
(NTE) [ 6 ]. The drawback of this approach is that the adjacency matrix at each time step is represented as the 
product of vectors, which is well-suited to highlight structure in communities but not adapted when the structure 
becomes more complex. 

Eollowing [ 20 ], we propose to use NME to find patterns in spectra of the collections of signals, obtained from the 
transformation of the temporal network. By analogy with music analysis, where an audio sample is decomposed 
into several audio samples, separating for instance voice from the instrumental part m, we would like to decompose 
the temporal network into temporal sub-networks, decomposing at each time step the global structure into several 
substructures. Eurthermore, audio spectra share similarities with graph spectra. We then propose to use the 
Itakura-Saito divergence as measure of dissimilarity, given by dis{x\y) = | ~ log y ~ allows taking advantage 
of the framework developed to reconstruct the signals from spectra as well as using efficient optimization algorithms 
with possible temporal regularization. 

As the input in our case is the temporal spectra 5, represented as a tensor of dimension C x F x T, a small 
adaptation has to be performed before applying NME. At each time instant t, the spectra St is represented as a 
vector Vt by successively adding end-to-end the columns of the matrix Sf. Eor alH G {0,... ,T — 1}, these vectors 
compose the columns of the matrix V, of dimension (FC) x T. The number of components K is set according to 
our expectations about the data, and the parameter 7 is strictly positive to ensure smoothness in the activation 
coefficients. 

5.3 Identification of components 

NME returns two matrices W and H : each column of W represents the kth (normalized) frequency pattern, while 
the kth column of gives the activation coefficients of the frequency pattern k at each time step. Erom a 
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Components Components Components 


(a) Frequency patterns, obtained after reshaping the columns of W into matrices. 



(b) Activation coefficients, corresponding to the rows of the matrix H, normalized by the maximal value of H 

Figure 6 : Results of the nonnegative matrix factorization for the Toy Temporal Network, using K = 3 and 7 = 6 . 


component-frequency map can be built by reshaping the vector into a matrix. To highlight how these structures 
are arranged in the temporal network, each component is transformed into temporal network: As described in 
m, using NMF with the Itakura-Saito divergence provides means of reconstruction of the collection of signals 
corresponding to each component. For each component /c, a temporal spectrum G is obtained using 

a Wiener filtering such that its elements read as: 


(k,t) _ '^{cf)khkt 

T^i=i'^{cf)ihit 

leading to a conservative decomposition of the tensor S: 

k=l 


( 10 ) 


( 11 ) 


The temporal spectrum of the component k is then a fraction of the original temporal spectrum. From 
an inverse Fourier transformation is performed, leading to a collection of signals for each component k denoted 
by G Finally, the adjacency tensor describing the temporal network corresponding to the 

component k is obtained by using the inverse transformation 7 ^ described in Section 


5.4 Application to the Toy Temporal Network (TTN) 

NMF is applied to the Toy Temporal Network defined in Section The matrix V G M^qsoxso decomposed 
using K = 3 and 7 = 6 , which is the expected number of structures. Figure 6 b shows the activation coefficients 
of each component. We notice that the activation coefficients are consistent with the division of time introduced 
in the generation of the temporal network: All components have distinct levels of activation, corresponding to the 
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(c) Component 3 
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(b) Component 2 


0 20 40 60 80 

(a) Component 1 


Figure 7: Adjacency tensor aggregated over time: for each component /c, 


four different structures used. The component 1 is active in periods 1, 2 and 4 with an almost constant level, the 
component 2 is mainly active in period 3, as well as in periods 1 and 4, and finally, the component 3 is active in 
period 1 and in period 2. 

Looking at the corresponding frequency patterns in Figure and more precisely their connections with the 
network structures observed in Section [231 confirms the good adequacy with the expected results: The component 
1 looks like a structure in communities, the component 2 resembles a /c-regular structure and the component 3 
exhibits random structure. The structure in period 1 is then a mixture between a random structure and a structure 
in communities (as it happens, one single community), in period 2 only the structure in communities is present, in 
period 3 a regular structure described by component 2 and finally, the fourth period is a mixture between structure 
in communities and regular structure. Random structure is present in period 1 and in period 2 inside communities. 

Figure shows, for each component /c, the aggregated adjacency matrix over time, obtained after the back 
transformation of the spectra into a temporal network. The sum is weighted by the activation coefficients given by 
the matrix iT, in order to highlight visually the most significant patterns: 

T 

/ifet (12) 

t=i 

This representation partly confirms the connections between spectra and structures as described above. We can 
notice that component 1 displays the three communities present in period 2, as well as communities corresponding 
to a regular structure in period 4. The actual communities in period 4 are caught in the component 2, which does 
not correspond to a regular structure, even if the diagonal, characterizing the /c-regular lattice, is clearly dominant. 
Finally, the component 3 looks like a random matrix, as no structure is visible, at least through this representation. 

The decomposition of temporal networks using NMF enables to retrieve from the spectra the different structures 
composing the temporal network, and to detect when these structures are present, either alone or jointly. We now 
propose to apply this framework to study a real-world temporal network describing contacts between individuals 
in a primary school. 

6 Temporal Network of Social Interactions Between Children in a Pri¬ 
mary School 

6.1 Description of the temporal network 

The decomposition is applied to a real-world temporal network, describing social interactions between children in a 
primary school during two days in October 2009. During a 20-second interval, an edge exists between two individuals 
if a contact is recorded, measured by wearable RFID (Radio Frequency IDentification) sensors [32]. For our study, 
the data set is described by a temporal network, representing for each time step the aggregated contacts between 
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Figure 8: Descriptors of the Primary School Temporal Network over time. The shaded regions correspond to 
class periods, while white regions correspond to breaks and lunch, according to the information given in [32]. No 
significant information are provided by the usual network-based descriptors 
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Figure 9: Marginals of the energies of temporal spectra.The energies of the low frequencies and of the first compo¬ 
nents are not equally distributed over time, indicating changes in the global structure of the temporal network. 


individuals for 10 minutes. We restrained the analysis to the first day: 226 children and 10 teachers participated in 
the experiment, separated in five grades (from 1st grade to 5th grade), themselves separated in two classes. 


show that the number of edges in the temporal network, as well as the number of isolated vertices, is not constant 
over time, reflecting the real-world nature of data: during the lunch break, some children leave the school to have 
lunch at home. The network-based descriptors (Figures and Figures [8^ do not provide significant insights about 
the structure of the temporal network. 

6.2 Spectral analysis 

Figure shows the marginals of the energies of temporal spectra obtained from the transformation of the primary 
school temporal network. The energies of the low frequencies and of the first components are not equally distributed 
over time: this indicates changes in the global structure of the temporal network, that occurs in break periods as 
well as during lunch. We can then divide the day into two main periods: the period where the children have class 
(shaded regions) and the breaks and lunch (white regions). 

Comparing the Primary School Temporal Network with a network with communities using correlations between 
spectra shows that the former temporal network is correlated with the latter network involving a large number of 
communities (between 9 and 15) during class periods, a smaller one (between 3 and 6) during breaks and lunch 
periods (see Figure pT]). This is consistent with the spatiotemporal trajectories given in [32], showing the location of 
the classes over time: during class periods, the classes are separated into different classrooms, while during breaks 
and lunch, the classes mix, yet in two distinct groups. 

6.3 Decomposition into sub-networks 

In the light of the above study of the temporal spectra, we can go further and decompose the Primary School 
Temporal Network into two temporal sub-networks. Furthermore, this allows for a quicker interpretation of the 
obtained components, as well as for an evidence of its ability the demonstration of the to extract the significant 
components. 

Figure [TT] shows the results of the NMF on the Primary School Temporal Network, using K = 2 and 7 = 5. 
The school day is divided into three specific periods, according to the activation coefficients (Figure [TT^. The first 
period occurs during class hours, where only the component 1 is mainly active. The second period groups together 
the breaks, during which the components 1 and 2 are significantly present. Finally, the period 3 concerns the lunch 
break, where the component 2 is dominant. This cutting is consistent with the ground truth, as described in [32]: 


Figure [^ shows the network-based descriptors used to characterize the temporal network. Figures [8a| and |8b 
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Figure 10: Correlation between the temporal spectra at each time step and the spectra of a network with a network 
with communities. Each line represent the number of communities, averaged over 20 repetitions. The temporal 
network is correlated with structure with a large number of communities during class periods, and with structure 
with a small number of communities during breaks and lunch. 
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(a) Frequency patterns, obtained after reshaping the columns of W into matrices. 



(b) Activation coefficients, corresponding to the rows of the matrix Ff, normalized by the maximal value of H 


Figure 11: Results of the nonnegative matrix factorization for the Primary School Temporal Network, using K = 2 
and 7 = 5. The activation coefficients are consistent with the schedule of the children in the primary school, as 
detailed in [3^ . 
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only two or three classes have breaks at the same time, while the other ones stay in their respective classrooms. As 
for lunches, they are taken in two consecutive turns, preserving the structure in classes, with nevertheless a weaker 
intensity. 

Figure \T2\ shows different representations of the temporal networks reconstructed from the components. The 
left figure shows the aggregated adjacency matrix over time, as described in Section [5^ The vertices are ordered 
according to the classes of the children, from the youngest to the oldest. The middle figure shows the aggregated 
network, using the layout provided in [32], after thresholding of edges according to their weights. The color of dots 
indicates the grade of the children, while black dots represent the teachers. Finally, the right figure shows the contact 
matrix between classes, obtained by counting the number of edges inside and between the classes. A logarithmic 
scale is used to enhance the visualization. Figure [T^ shows the original temporal network, while Figures [T^ and 
|12c| show respectively the component 1 and the component 2. We can easily observe that the component 1 describes 
the structures in classes, with higher density of edges inside classes than between classes. Conversely, the component 
2 highlights a less structured network pattern, which looks like two communities, separating the youngest classes 
from the oldest. Those observations are consistent with the description of lunches mentioned above. 

These results highlight the interest of decomposing a temporal network into several sub-networks, which can be 
studied independently of one another. Without prior knowledge, the different periods of activity in the primary 
school are displayed, and can then guide the analysis of the system by restricting the analysis over several time 
intervals. 


7 Conclusion 

We have proposed a novel method to track the structure of temporal networks over time using the duality between 
graphs and signals, as well as classical signal processing techniques, such as spectral analysis and nonnegative matrix 
factorization. At each time, the temporal network is represented by a graph which is transformed into a collection of 
signals. NMF is used to extract patterns in the energies of spectra of these signals; these patterns, whose activation 
coefficients vary over time, represent a specific structure of the underlying network. The effectiveness of the method 
has been demonstrated on a toy example, containing three types of structures as well as on a real-world network 
describing temporal contacts between children in a primary school. 

These results provide insights in the characterization of temporal networks, but also call for further studies: In 
particular, it would be interesting to have a deeper look on how the reconstructed temporal networks are embedded 
in the original temporal network. Furthermore, it would be worth considering an iterative analysis, by studying the 
reconstructed temporal network themselves, using the same process. This could lead to a spatiotemporal multiscale 
analysis, permitting an increasingly precise track of the structure of the temporal network. 
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Figure 12: (Left) Aggregated adjacency matrix over time, weighted by the coefficients of H. (Middle) Network 
representation using the layout provided in [32], after thresholding of edges according to their weights. The color of 
dots indicates the grade of the children, while black dots represent the teachers. (Right) Grayscale-coded contact 
matrix between classes: each entry gives the number of contacts inside and between the classes. A logarithmic scale 
is used to enhance the visualization. The component 1 represents the structures in classes, while the component 2 
describes the structure during the breaks and lunch. 
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